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We study a system of interacting self-propelled particles whose walking velocity depends on the 
stage of the locomotion cycle. The model introduces a phase equation in the optimal velocity 
model for vehicular traffic. We find that the system exhibits novel types of flow: synchronized 
free flow, phase-anchoring free flow, orderly jam flow, and disordered jam flow. The first two 
flows are characterized by synchronization of the phase, while the others do not have the global 
synchronization. Among these, the disordered jam flow is very complex, although the underlying 
model is simple. This phenomenon implies that the crowd behavior of moving particles can be 
destabilized by coupling their velocity to the phase of their motion. We also focus on “phase¬ 
anchoring” phenomena. They strongly affect particle flow in the system, especially when the density 
of particles is high. 


I. INTRODUCTION 

The movements of all animate things, such as walking 
of humans, flapping of birds, and swimming of fish, are 
restricted by their locomotion jj'. Since the collective 
motion of self-propelled particles (SPP) began to attract 
interest , many studies have been conducted in vari¬ 
ous contexts @, including pedestrian crowds m , bird 
flocks [l5|, and insect swarms EM3- However, the ef¬ 
fect of locomotion on macroscopic behavior remains an 
open question. In fact, recent experiments in the field 
of pedestrian dynamics suggest that the oscillatory mo¬ 
tion of pedestrian walkin g m ay have a significant effect 
on overall dynamics [H| |lj]. Jelic et al. studied the 
collective dynamics of pedestrians walking in a line, fo¬ 
cusing on bilateral oscillations during walking [l3|. From 
fluctuating trajectories, they obtained the phase of loco¬ 
motion and investigated the interaction between neigh¬ 
boring pedestrians. Furthermore, Yanagisawa et al. [14j 
reported that synchronization of the locomotion phase 
among pedestrians may increase the total flux. Here the 
“phase” is a mapping from each stage of locomotion to 
a real number if> £ K/27tZ. Unlike the situation in the 
vehicular traffic, this phase is strongly related to the ve¬ 
locity of walkers. 
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FIG. 1. Space-time diagrams (top) and velocity-time dia¬ 
grams (bottom) for orderly jam flow (left, p = 2.0, K = 
5.0, A = 0.05) and disordered jam flow (right, p = 3.0, K = 
1.0, A = 0.05). Bold lines indicate the trajectory of the 50th 
walker. The small regular waves in orderly jam flow corre¬ 
spond to phase anchoring. 

In this paper, we propose a simple model that describes 
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particle-following behavior, phase-velocity coupling, and 
phase-phase coupling. Particle-following behavior is im¬ 
plemented by the optimal velocity (OV) model 0 for ve¬ 
hicular traffic, which is known to emerge the phase tran¬ 
sition from free flow to jam flow that is highly tractable in 
mathematical analysis. Phase-phase coupling is consid¬ 
ered within the framework of the Kuramoto model [l9l ] , a 
paradigmatic model describing synchronization phenom¬ 
ena in nature. In the present study we assume that the 
phases of successive particles tend to synchronize, which 
is experimentally suggested when pedestrians walk in a 
dense crowd [Hj. Since each walker follows its predeces¬ 
sor and ignores its successors, the interaction is unidi¬ 
rectional and local. In contrast, other studies have been 
mainly devoted to globally coupled oscillators or locally 
but bidirectionally coupled oscillators [2(J[2l|. Note that 
we do not consider detailed modeling of actual animals; 
instead, we concentrate on ideal particles called oscilla¬ 
tory walkers (OWs). Although here we assume the in¬ 
teraction between particles is through their headway dis¬ 
tances, it could depend also on their (relative) velocities, 
which is beyond scope of this paper. 

In spite of its simplicity, the OW model presents a rich 
behavior, including a novel type of jam [Fig. [2(b)]. We 
report its fundamental characteristic and give insight into 
collective locomotion. 

The rest of the paper is organized as follows: The next 
section is devoted to define the model. In Sec. uni we 
summarize the dynamics of the model. To understand 
its destabilization phenomena, theoretical analyses are 
performed in Sec. IIVI Finally we summarize these results 
in the final section. 


II. MODEL 

Consider N particles, labeled n = 1, ■ • •, N from left, 
walking to right on a one-dimensional ring. We impose 
periodic boundary conditions [particle (N + 1) particle 
is identical to particle 1]. The velocity of the nth walker 
(x n ) is determined by the following equation: 

x n = a{V(Ax n ) + A(cosc/) n + 1) - x n }, (1) 

where V(Ax n ) is the optimal velocity determined for 
headway distance ( Ax n = a; n +i — x n ) [l8|. Each walker 
adapts its velocity to this optimal velocity with an adap¬ 
tation intensity (the reciprocal of the reaction time) a. 
The effect of locomotion is represented by the oscillation 
term A(cos <f> n + 1). The quantity A is the magnitude of 
the fluctuations in the target velocity, which is set to a 
small value. The fluctuation phase (j> n is determined by 

4> n = ui(Ax n ) + A'sin A(j) n . (2) 

Here we assume that the ideal angular velocity u>(Ax n ) is 
defined by a function similar to the ideal walking velocity 
V(Ax n ), i.e., using a normalized optimal velocity (NOV) 
function U(Ax n ), we can express uj(A x n ) = QmU(A x n ) 


and V(Ax n ) = VmU(A x n ) . We further impose the con¬ 
dition that when the headway distance is large enough, 
particles walk with a steady angular velocity HmC^oo and 
velocity VmUqo , and for small Ax n , U(Ax n ) ~ 0 to 
stop. To capture these assumptions, the NOV function 
is set to be a differentiable and monotonically increas¬ 
ing function that is asymptotic to 0 as Ax n —>• 0 and to 
a positive constant as Ax n —» oo. We use the follow¬ 
ing form of the NOV function when actual calculations 
are needed: U(Ax n ) = ci[tanh(c 2 Ax ra — C 3 ) + tanhcs] 
with the positive parameters ci,C 2 and C 3 . U(Ax n ) is a 
monotonic increasing function such that 17(0) = 0 and 
U{ 00 ) = Ci [1 +tanhc 3 ]. The quantity C 3 /C 2 corresponds 
to its inflection point, and c\ and C 3 determines its height 
given by Ci[l + tanhc 3 ]. Here we set c\ = 0.5, C 2 = 5, and 
c 3 = 2.5. 

The second term on the rhs of Eq. @ represents the 
synchronous interaction of successive walkers, whose in¬ 
tensity is AT > 0. When the phase difference A cj) n = 
4> n +i — 4> n is zero and the walkers are uniformly dis¬ 
tributed, they walk with a steady velocity and rhythm 
(synchronized free flow). In this paper, we restrict the 
model parameters to a = 3, Vm = 1 ,Hm = 1 and 
N = 100. We vary the particle density p = N/L by 
changing the system length L. 


III. DYNAMICS OF THE OW MODEL 

The OW model has four types of flow: synchronized 
free flow (SFF), phase-anchoring free flow (PFF), orderly 
jam flow (OJF), and disordered jam flow (DJF). “Phase 
anchoring” occurs when phase differences A <f> n are fixed 
at a common value, A^o G (—When the effects 
of distance between walkers are negligible, the model re¬ 
duces to the Kuramoto model with local interactions. 
Then, the states A<p n =const., in which all angular ve¬ 
locities are identical, are stable over a certain range of 
A(j) 0 . This property is inherited by the OW model; how¬ 
ever, as can be easily verified, the states A <j> n =const.^ 0 
cannot be realized exactly because of the effects of spa¬ 
tial structure. Actually, each phase difference fluctuates 
slightly around Ac/)q as depicted in Figs. [2(c) and (2d). 

Before going into details, we briefly summarize the dy¬ 
namics in the system. First, similar to vehicular traffic 
flow, the model has phase transitions between free flow 
and (orderly) jam flow; these are caused by increases 
in the density of walkers. When the density is small 
{p < p C r,i), free flow is stable and a jam cannot develop. 
At middle densities p cr> i < p < p C r, 2 , free flow is no 
longer stable, and small perturbations grow into a jam; 
this is called OJF [Figs. [2a) and [2(e)]. At high densi¬ 
ties Per ,2 < Pi the system regains stability for free flow. 
In the free flow regime, the OW model has two types of 
flow: SFF [Figs. [2(a) and [2(b)] and PFF [Figs. [2(c) 
and [2d)]. These flows are locally stable, and all tran¬ 
sient flows are asymptotic to either SFF or PFF in the 
low-density regime, depending on the initial condition. 


3 




■H 

0 




-Ji/2 

(a) SI 

p=\.0, K=5.0 

H F (low density) n 

1 

•H 


0 


7T/2 


(b) si 

p=3.0, A=5.0 

H F (high density) " 



(d) PFF (high density) 


1 



n/2 

~n/2 

(e) O. 


p-2.0, K= 5.0 

IF 

1 

-H 


7l/2 

-n/2 

/j-3.0. A-5.0 


(f) DJF 


FIG. 2. Snapshots of velocity and phase difference for A = 0.05. On the abscissae, n represents particle number. Panels (a) 
and (b) are for synchronized free flow (SFF), panels (c) and (d) are for phase-anchoring free flow (PFF), panel (e) is for orderly 
jam flow (OJF), and panel (f) is for disordered jam flow (DJF). Wavy shapes of velocity in (c), (d), and (e) stem from the 
constant phase difference. In PFF (c) and (d), the phase differences are almost constant, but slight fluctuations occur, as shown 
in the insets. For a = 3.0 we set, p = 1.0, 2.0 and 3.0 correspond to low-, middle-, and high-density regimes, respectively. As 
explained in the text, the middle-density regime has OJF only. 


TABLE I. Classification of flows. Each tick indicates the ex¬ 
istence of the flow in the density regime. 



Density 

Low 

Middle 

High 

Free flow 

Synchronized free flow (SFF) 

/ 


/ 

Phase-anchoring free flow (PFF) 

/ 


/ 

Jam flow 

Orderly jam flow (OJF) 


/ 


Disordered jam flow (DJF) 



/ 


In contrast, in the high-density regime, SFF is always 
locally stable, while PFF can be destabilized for some 
parameter regimes, resulting in another type of jam. 

In OJF of the middle-density regime, separation of 
phase differences occurs. As shown in Fig. [2e), the 
phase differences become almost constant in each region. 
This phase anchoring is produced by the spatial struc¬ 
ture of OJF, and it also depends on K and the initial 
conditions. 

In addition to these states, another nontrivial phe¬ 
nomenon emerges in the high-density regime, namely dis¬ 
ordered jam flow (DJF). This is complex jam flow in¬ 
duced by strong coupling between <j) n and x n . As shown 
in Figs, [ijb) and 0)0 , no steady pattern occurs, and 
small jams propagate with different velocities and even 
disappear and reappear. The four types of flow correlate 
with density as in Table. [TJ 

Now, we investigate the system in detail. 


A. Phase-anchoring phenomena 


We first focus on phase-anchoring phenomena. In this 
paper, we use the following definition for PFF: (i) Each 
Acj) n and Ax n periodically fluctuates around a certain 
phase difference Afo (const.) and headway distance 
Axq = respectively, (ii) The ranges of these fluc¬ 
tuations are small (<C |A(/>o|). 


The effects of phase anchoring are significant, espe¬ 
cially when the density of walkers is high. In this regime, 
since each Ax n and ui{Ax n ) is small, the dependence of 
the system dynamics on phase becomes relatively large. 
One intriguing feature is the relationship between phase 
anchoring and the flux of walkers. We show the relation 
between flux and phase difference Atfio in Fig. [3] Flux 
is defined as the product of density and average velocity, 
averaged over time, J = {pi). For SFF, we can calculate 


the flux as J s = p 




A 


Simulation results indi¬ 


cate that for small phase delays and advances, the flux 
is decreased, but for large phase differences, the flux in¬ 
creases, especially when the phase advances (A cj> 0 < 0). 
In addition, as expected, the effect of phase difference is 
conspicuous when the intensity of the coupling parameter 
A is large. Although phase-anchoring can be observed in 
the low-density regime, there is little change in the flux 
because when the headway distance is large, fluctuations 
in phase anchoring do not influence the angular velocity 
through the term tu(Ax n ) ~ VLm- In PFF, each Ax n 
oscillates with a certain phase that differs from the walk¬ 
ing phase </>„. This difference distorts the trajectory of 
walkers from that in SFF (represented by a sine curve), 
leading to decreases and increases in the flux. (See also 
Appendix iBl) 
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FIG. 3. Relationship between flux and phase difference for the 
case of p = 4.0 and K = 5.0 (PFF in the high-density regime). 
Here the flux is normalized by the SFF value (A0o = 0). 



Ax\ 


FIG. 4. Headway-velocity trajectory for DJF (A = 0.05, p = 
3.0, and K = 1.0). The red broken line represents the trajec¬ 
tory of jam flow in the OV model (A = 0, p = 2.0). 


B. Disordered jam flow 

For small K and high density p, free flow allows only 
synchrony (A</>o = 0) and PFF becomes unstable, and a 
highly complex flow appears. This flow is phenomeno¬ 
logically distinct from OJF, which is already observed 
in the OV model. Instead, it originates from the strong 
coupling between x n and 4> n . In DJF, small jams locally 
appear and disappear intermittently, and no stationary 
state is reached [Fig. [TJb)]. Further, we can no longer 
observe regular patterns in x n , x n , or A<fr n [Fig. [ljb) and 
EJf)]. Figure |4] shows a complex limit cycle of DJF in the 
(Axi,i:i)-plane; the figure also shows the trajectory of 
the jam flow in the OV model. In the quasi-periodic tra¬ 
jectory of DJF, we can see distinctive patterns, especially 
at its edges. 


IV. THEORETICAL ANALYSES 

Next, we discuss the stability conditions for each flow. 
To understand the destabilization of steady free flow, 
(a;^, ), we investigate how the flow responds to small 

perturbations. For a small perturbation, x n = Xn ^ + 


= <j)n + S n , the system equations m and m lin¬ 
earize to 

i n /a = -i n + V M U'(Ax^n^Acn - A sin </4 0) <5„, (3) 
Sn = K cos A4 0 ) A 5 n + n M U'( Ax^)Ae n (4) 

with U'(Axn^) = When all walkers are 

uniformly synchronized, namely, Asi 0 * = L/N = 1/p, 
<; t>n ^ = fl M t/(l/p)t = (j> s t, Eqs. ([3]) and 0) simplify to 

in/a = -in + V M U'(l/p)Ae n - A sin (<j) s t)5 n , (5) 

5 n = KA6 n +n M U\l/p)Ae n . (6) 

Here we assume that <p s t is a fast variable compared to 
the growth rate of the perturbation, i.e., e, S, and their 
derivatives are constant for each period T = By 

integrating Eq. over this period, we can eliminate 
the vibrational term, leaving 

in/a = -i n + V M U' {1/p)Ae n , (7) 

which is the linearized OV model. Now, the stability 
condition is given by VmU'( 1/p) < | [22j]. When e n does 
not grow in Eq. 0 , equation 0 gives the condition on 
K as I< > 0. 

For PFF, the analysis is not straightforward; we start 
with Eq. ©. By ignoring the coupling between e n and 
6 n as well as the fluctuations in the phase difference (i.e., 
A (jin'* ~ A</>o), we obtain 

Sn = K COS ( A(j) 0 )ASn■ (8) 

From this equation, we can find the stability con¬ 
dition, | A</> 0 1 < Additionally, since the system 

is periodic (</>n+i = <f>i)i the condition NAtfig = 
2 ?n 7 r (m = ±1,±2, ■••) must be satisfied. To summa¬ 
rize, stable A(f>g can take the discrete values of A(f>g = 
0 , ±^r, ± 7 ^-, • • •, ± ^ where is the floor 

function giving the maximum integer not larger than Y. 
These conditions were verified by numerical simulations. 
On the other hand, the coupling behavior of e n and S n 
is determined by the final term in Eq. 0. If the time 

integral f^ sin ^ 0 (i.e., if the sine function is 
“distorted”), the system can destabilize. This distortion 
is also observed as a change of flux in the high-density 
regime (Fig. [3] See also approximate analyses in Ap¬ 
pendix [Bj). From these facts, we conclude that the high- 
density regime leads to coupling in the stability equation, 
which may result in DJF. 

In Fig. [5]we illustrate the phase diagrams of the system 
by simulations, using the expression for PFF (x n ,i, 4>n, l) 
(whose derivation is found in Appendix!© as initial con¬ 
ditions, to see whether the flow is stable or not: 

n Aa . , 

x n ,i = vot H- 1 - , sin W (t) 

p <j) 0 \/a 2 + <j% 
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FIG. 5. Phase diagram in the high-density regime. Initial 
conditions are given by the expression of PFF with A0o = 
(x n ,i, 0n,i)- If we take SFF as initial conditions (A0o = 
0), the flow is always stable. Green diamonds indicate the 
parameter region where PFF is stable. Orange squares and 
red circles correspond to the emergence of DJF, and the latter 
include collisions of walkers. The same diagrams for low- and 
middle-density regimes only show PFF and OJF, respectively. 


2Aa 2 V'(j;)sm{^) 

^o(° 2 + </>o) 

4>n, i = 0o* + A0o n 

22W(i)sin^ 


sin(^(t) - ip 0 


A 0 O , 


sin(<I>(t) + 


where vo = V(l/p) + A, 0o = w(l /p) + K sin A0o, <Po = 
tan _1 (0o/a), and 'P(i) = 0o t — V^o + A0o n. It is found 
that the system destabilizes for small K and large A, 
and for especially large values of A , particles even collide 
with each other in DJF. Thus we can conclude that the 
increase of A and K are connected to the destabilization 
and stabilization of the system, respectively. In addition, 
increase of particle density reduces the DJF regime, be¬ 
cause it decreases space for “over acceleration” of parti- 

I 


cles that is indispensable for DJF. 


V. CONCLUSIVE DISCUSSION 

In summary, we have proposed a simple model to show 
that the restriction coming from locomotion largely af¬ 
fects the collective behavior of SPPs. The model has var¬ 
ious types of flow, including a novel type of complex flow 
called disordered jam flow (DJF). Disordered jam flow 
appears when the particle density is high; this flow is 
caused by coupling between velocities and phases of par¬ 
ticles. We also find that particle flux is strongly affected 
by phase anchoring. On the other hand, some issues re¬ 
main to be clarified in future studies. These include, 
synchronization phenomena on inhomogeneous distribu¬ 
tions of intrinsic optimal velocity and angular velocity, 
quantitative analysis for phase anchoring and its effect 
on flux, etc. These points should be addressed in detail 
in future works. 

Although pedestrian locomotion in a crowd is not yet 
fully understood, a real pedestrian may adapt the crowd 
velocity and step length by considering the environment 
of other pedestrians (not only the nearest predecessor). 
Furthermore, the environment may affect a pedestrian 
through a strong psychological repulsive force 0- We 
believe that by including the effects of phase, models for 
pedestrians and other animals can be improved. More¬ 
over, the synchronization of pedestrians through a bridge 
23, ;24] might be an interesting problem. It is also at¬ 
tractive to extend this model to two dimensions using 
the two-dimensional OV model [25|, which generalization 
would require additional rules of phase synchronization. 
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Appendix A: Approximate description of PFF 


In this section we give equations of a particle trajectory in PFF approximately. The trajectory satisfies the following 
system equations: 

x n = a{V(Ax n ) + A( cos0 n + 1) - x n }, (Al) 

0 n = uj(Ax n ) + K sin A 0 n . (A2) 

First, we assume that the deviation of each Ax n is small compared to its absolute value, and that the phase difference 
A 4> n is constant (A0o). In this assumption, Ax n = i,0„ = w(^) + /vsinA0o = 0o (const.). Then, the system 
equations m and (IA2I) can be easily integrated: 


4>n, o — 4>o t + A<f>o n, 

. . Aa ■ . n 

Xn,o{t) = - , sin (00* - +o + Acpon) + v 0 t + -, 

p 


<Wa 2 + 0o 


(A3) 

(A4) 
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FIG. 6. Relationship between Aa and Apo for the cases of A = 0.05. 


x n ,o{t) 


Aa 

\j a 2 +P o 


COS ((pot 


<Po + A tpon) + Vo- 


(A5) 


Here we put vq = V(l) + A and ipo = tan 1 (^-). These equations give exact trajectories of particles when A (p 0 = 0 
(SFF). However these expressions are not enough for understanding the dynamics of PFF in high-density cases. We 
next give a more accurate approximation based on these equations; the deviation of particle distance Ax n is taken 
into consideration, while the deviation of phase-difference is still ignored. Using the previous result we approximately 
give Ax n = x n+ % t o — x Ui o in advance: 


Ax n -~ 


Aa 


sin ((pot — (p 0 + A(p 0 n + Acp 0 ) - sin (p 0 t — ipo + A(p 0 n) 


Po\J a 2 + Po 

2 Aa sin Aia / ■ A <p 0 \ 

: cos ( (pot -ipo + A(p 0 n + —- ) 


a 2 + Pi 


= A a cos ( (p 0 t — ipo + A(p 0 n + 


A(p 0 


where Ha is the amplitude of oscillation of Ax n . Then, the equation for (p. 

(pn ,i = w(l/p) + K sin A^o + w'(1/p)Ha cos ((pot - (po + A(p 0 n + 


is now rewritten as 

A(po \ 


2 )' 


(A6) 

(A7) 

(A8) 


(A9) 


which can be solved as: 


Pn, 1 — 

Here we assumed Ax n — A <^C 1 

*^n, 1 — 


Dn, 0 


+ 


2Aau'(l/p) sin ^ 


p 2 o\/a 2 + (y 0 


sin p 0 t — <Po + A p 0 n + 


A 0 O ^ 


. In the same manner, we can derive a more accurate expression for Ax r 


, 2Ha 2 W(l/p)sin(^) . ^ o , A ,A0 O , 

x n ,o H---xr-sm (p 0 t - 2tp 0 + Ap 0 n + 

Po( a + Pi) 2 


(A10) 


(All) 


Note that the effect of these correction terms increases as |<fo| (or A^o) becomes small (large), namely, the first 
approximations presented in this section do not agree with actual trajectories. As shown in Fig. [6l the approximation 
is valid for large K and small positive A p 0 . 

We improved the approximation using (x U: o, pn, o) to obtain (x Ut i, p n ,i)- In the same procedure, further improvement 
can be made by calculating the series of trigonometric functions (x n ,k+i, pn,k+i) from (x n ,k, p n ,k) if their amplitude is 
small enough to allow the Taylor expansion of the optimal-velocity functions, while it is far more difficult to describe 
the system where this assumption is not valid. 

In the main text, we used x„ : i and p n ,\ as initial conditions to plot the phase diagrams. 
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(b) *0« 


FIG. 7. Relationship between (f> o and A<f>o for the cases of A = 0.05. 


Appendix B: Distortion of trajectory and its effect on particle flux 


Average velocity of a particle in a steady flow is given by integrating Eq. (I A 111 over a period of one “step”, T: 

j: j x n dt=Y j V(A x n )dt + A+Y J A cos <j> n dt. (Bl) 

By using the first-order approximation (x n fi, <j) n , o) and assuming that the amplitudes of Ax n and Aa are small enough, 
we can estimate the average velocity as 


Y J x n dt~v o + y J Acos^pdf. 


(B2) 


Therefore the particle flux in PFF is affected by the fluctuation of 0 n from 4>n ,o , by which the final integral in the 
above equation remains nonzero. Now the degree of the distortion is characterized by Eq. (IA9I) . For example, we 
see the difference between the functions cos 0 nj o and cos</> n> i when the focal phase difference Apo is a small positive 
number that leads to small 4>q. The former function is at its maximum at pot + Apon = 2rmr, at which the phase 
velocity of the modified approximation p ni o is faster than <f>: 


Z,1 \ < f ) n,0 = ^rn'K — 00 A 


2 Aau/(l/p)sin^a /A p 0 


■ COS 


00 v a 1 2 3 4 5 + 0i 


- To > 00 - 


(B3) 


The contribution of the phase difference, — tpo, is depicted in Fig. [ 8 ] In the same manner, we can find that when 
( t>o = (2m + l) 7 r, is slower than cf>. Consequently, the positive part of the cosine function (that corresponds to 
the approximation 0 nj o) shrinks while its negative part is amplified (see Fig. [9]), resulting in a smaller value of the 
integral, and the particle flux decreases as shown in Fig. 3 in the main text. Thus the decrease in the flux around 
A<po = 0 can be explained in terms of the interaction between oscillations. However, this theory cannot account for 
the cases of J/J s > 1, for which we have to take large deformations of trajectories into account. 
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